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Evidence from the lattice suggests that formation of a flux tube between a qq pair in the QCD 
vacuum leads to quark confinement. For large separations between the quarks, it is conjectured that 
the flux tube has a behavior similar to an oscillating bosonic string, supported by lattice data for the 
groundstate qq potential. We measure the excited states of the flux tube in 3d SU(2) gauge theory 
with three different couplings inside the scaling region. We compare our results to predictions of 
effective string theories. 



I. INTRODUCTION 

Simulations over the last few years have accumulated 
strong evidence that gluonic dynamics indeed leads to 
the formation of a flux tube between test quark and an- 
tiquark (qq) in the vacuum of Yang-Mills theory 0, Q. 
This implies a linearly rising potential between quark and 
antiquark in the QCD vacuum and thus leads to quark 
confinement. At large qq separations, this flux tube is 
expected to behave like a string. Open bosonic string 
descriptions of the dynamics of this flux tube have been 
attempted for a long time J3] . Using the Nambu-Goto 
(NG) action, first Alvarez [J] (in the limit d — » oo) and 
later Amis [5| obtained the energy states of the flux tube 
as 



E n (R) =o-R 
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(d-2) (1) 



where a is the string tension, R the quark-antiquark sep- 
aration and d the number of space-time dimensions. A 
closed string description was proposed by Polchinski and 
Strominger (PS) p] where they suggested how effective 
string theory with vanishing conformal anomaly could 
be formulated in arbitrary dimensions. The spectrum 
of the string using the PS prescription has been com- 
puted in @,H, Q. To order R~ 3 it has been found that 
the spectrum is universal (depends only on the number 
of space-time dimensions and the string tension), and to 
this order it coincides with the NG spectrum. For a calcu- 
lation of the closed string spectrum on the lattice see [l(| • 
Using symmetry properties of the Polyakov loop correla- 
tion function, Luscher and Weisz (LW) showed that any 
effective string picture posseses an open-closed duality 
property They also showed that demanding this 

duality constrains the possible string spectra and in par- 
ticular forbids 1/R 2 terms in the effective string action. 
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In the LW formulation too, the spectrum is consistent 
with the NG spectrum. In 3-dimensions it is exactly the 
same and in 4-dimensions one undetermined parameter 
remains which can of course assume the NG value. The 
NG partition function itself respects open-closed duality. 
Predictions from these theories can be tested by com- 
paring them to results coming from simulations of pure 
Yang-Mills theories. 

In this Letter we will present a new scheme for mea- 
suring observables in simulations of pure Yang-Mills the- 
ories. Using this scheme we will accurately measure ob- 
servables in 3-dimensional SU(2) lattice gauge theory 
and try to extract the excitation spectrum of the flux 
tube. For earlier studies of the QCD string spectra using 
different schemes in both 3- and 4-dimcnsions for several 
gauge theories, see [T3| and the references therein. 



II. PRELIMINARIES 

To probe the properties of the flux tube formed be- 
tween test quark and antiquark, the two observables in 
pure Yang-Mills theories are the Polyakov loop correla- 
tion function and spatio-temporal Wilson loops. While 
Polyakov loop correlation functions project very strongly 
onto the ground state, Wilson loops are much more sensi- 
tive to the excited states of the flux tube and in addition 
allows one to couple more strongly to a particular state 
while suppressing the others, by appropriate choice of the 
spatial parts of the loop. It is therefore the observable of 
choice if one wants to study the excitation spectrum of 
the flux tube. 



A. Extraction of the excited states 

In 3-dimensions, the states of the oscillating string can 
be classified by charge-conjugation and parity properties 
(C, P) [H, Choosing the spatial parts of the Wilson 
loops as shown in figure[TJ the (C, P) projectors are given 
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by the superpositions 

S++ = Si + S 2 + §3 + §4 
S+- = Si +§ 2 - § 3 - §4 

(2) 

S~ = Si - s 2 - s 3 + s 4 

S _+ = -Si + S 2 - §3 + §4 

where +/— stands for even or odd states under C and 
P. Although in each of the (C, P) channels, there are an 
infinite number of states, we are going to look only at the 
ground states and will label these channels as {0, 1, 2, 3} 



respectively. The Wilson loop projecting onto a channel 
n has the spectral representation 



W n {R,T) = Y.W) e ~ K(R)T - ( 3 ) 

i=0 

where E % are the energies of the states in a given chan- 
nel. Using Wilson loops with different temporal extents, 
one obtains the energies and energy differences at leading 
order as 
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J-b - J- a v 1 

AE nm (R) + -1— a n (R) e~ s ^ T « (l - e"^ ^ 
J-b — J-a v 



(4) 
(5) 



where T a < T b . In the following, the values obtained 
from the LHS of eqns.JH) and (JSJ) will be called naive 
values and be denoted by E n (R) and AE nm (R). The 
quantities E n (R) and AE nm (R) will be called corrected 
values. They, along with a n and S n , are obtained as 
fit parameters. The fits are done by taking all possible 
combinations of T a and T b and are discussed in more 
detail in Section D. 
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FIG. 1: Top: Basic set of operators used to construct (C, P) 
channels. Bottom: Set of operators with stronger coupling 
to the second excited state. 



B. Noise reduction techniques 

String like behavior of the flux tube is expected to oc- 
cur at large qq separations. The Wilson loops we measure 
must therefore extend to large enough R values. Also 
we see from equations (|4]) and ([5]) that to reduce the 
contaminations due to other states we must either go 
to large values of T or tune a to small values by ap- 
propriate choice of the basis states. The latter method 
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has been followed in where asymmetric lattices with 
small temporal extents were used. The former method 
requires accurate measurements of expectation values of 
large Wilson loops, made possible recently by the multi- 
level algorithm, proposed in [161 ]. This method has been 
used in [It], EH to accurately measure the ground state 
properties of the flux tube. Here we will try to combine 
both methods by using a slight variant of the multilevel 
algorithm. 

In the multilevel algorithm, intermediate averages are 
computed for the temporal links by updating certain sub- 
lattices with the sources on the spatial links of the bound- 
aries of these sub-lattices. The boundaries are held fixed 
during the sub-lattice updates. We now put the sources 
on space-like surfaces in the middle of a sub-lattice of 
thickness 2a where a is the lattice spacing. Two tempo- 
ral links are attached to the two ends of the source which 
terminates on the space-like surfaces that are fixed during 
the sub-lattice updates (see fig© • The advantage is that 
the sub-lattice updates reduces fluctuations of the sources 
as well as fluctuations of the temporal links. Moreover, in 
addition to using multihit on the temporal links, we can 
now use multihit on the spatial links as well. For further 
details see I13L Il9j [ and references therein. Adopting the 
notation in we now define the expectation value of 
the Wilson loop by 

(W(T)) = <{L(0)} Q7 {T( a ) ■ ■ • T((T-l)a)} QW {L(T)*}^), 

(6) 

where L(t) is the operator of fig(2]and { } indicates sub- 
lattice averaged quantities. The only difference from [l6| 
is that now we have {L} instead of L. A similar method 
was used in Q to measure the static potential and ob- 
serve string breaking in the adjoint representation. 
In most parts of the calculations, the operators {S/ 1 } 
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FIG. 2: The source of the Wilson loop in the modified algo- 
rithm. The product of the directed lines defines the source. 
The thick lines are held fixed during sub-lattice updates. 

(figure[T]top; see also [HLGjl]) are good enough for a reli- 
able signal. However, for the second excited state beyond 
/3 = 7.5 the error reduction was not sufficient and we had 
to use another set of operators {§f }, (figure Q] bottom) 
with a stronger coupling to that particular excited state. 
We will henceforth refer to the two sets as set A and set 
B, respectively. To check the error reduction with the 
new set of operators we performed 460 measurements at 
(3 = 7.5 with spatial extents R = 15 — 20 using the simu- 
lation parameters shown in table HI The resulting relative 
errors of Ei are also tabulated there. The values for the 
naive energies using both sets are plotted in figure [3J For 
comparison we also plot £2, which was calculated from 
E2 , in the same region. It is clearly seen that not only do 
the naive values of set B already coincide with the cor- 
rected values of set A, they also have much smaller error 
bars. Set B therefore significantly improves the signal 
for the state S~~. 

There are several parameters for the algorithm, that 
have to be tuned properly in order to achieve the maximal 
error reduction. For the temporal links we follow the 
same procedure as outlined in In addition we now 
have another parameter which is the number of updates 
(N s ) for the sub-lattices containing the spatial operators. 
We found it beneficial for the excited states to have quite 
a large number of such sub-lattice updates. Since there 
are only two such sub-lattices for every loop this does not 
increase the cost of the simulation very much. Table |TT] 
contains the resulting optimized run parameters and we 
refer to [l3| for details of the optimization. 

C. Simulation parameters and lattice scales 

Our simulations were done with three different cou- 
plings, (3 = 7.5, 10.0 and 12.5 (chosen so as to lie in 
the scaling region) in 3-dimensional SU{2) lattice gauge 
theory, using usual heatbath sweeps (2lT |. combined with 
three overrelaxation sweeps. The scale was set by the 
Sommcr parameter ro |22j , which has been computed for 
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TABLE I: Left: Parameters of the testruns to compare the 
operators. Right: Relative error of E2 in % for operator sets 
A and B. 




FIG. 3: Comparison of the naive energy E2 from set A 
and B. 

these (3 values with high accuracy in [l8(. For each beta 
value we calculated Wilson loops with four different tem- 
poral extents and used cubic lattices. 

Our simulation parameters, lattice scales, along with 
the parameters for the multilevel algorithm arc also tab- 
ulated in table ITU 

D. Error analysis and control of the fits 

For estimating the error of all the energy values and 
differences, we used the usual binned jackknife method, 
with 44, 24, 40 and 40 bins for the lattices L A , L% L% 
and L , respectively. We also checked that the errors did 
not vary by more than a few percent with bin size. 

For the energy difference AE 2 o at (3 — 10.0 and 12.5 
we used the results obtained from set B for E2 while the 
values for Eq were obtained from set A. Since these two 
values come from independent simulations, we added the 
individual errors in quadrature to obtain error estimates 
for Ai?2o- 

The remaining issue is the control of the fits (0| and 
([5]), and this was done in two ways. For the energies E n , 
we expect a to be smaller then the ratio of the degen- 
eracies [13] of the energy states considered and S should 
be of the order of the energy gap to the next level in the 
channel. Similar conclusions hold for the parameters of 
the energy differences [l3|, [TiJ as well. If the resulting 
fit parameters were far of from these expectations we did 
not trust the fits. As a second check we plot the expected 
corrections 

A = ^-L- a n (R) e- s ^ T a ^ _ e - Sn (R) (t 6 -t^ 

obtained with averaged parameters a n and 5 n against the 
differences A = E n — E n , for all possible combinations 
of T a and T^. We show an example for these fits (AEiq, 
(3 = 7.5, R = 11 and 15) in figure |U If for some of the 
values there was a big discrepancy to the expectation 
A = A the fit was also regarded as unreliable (for more 
details and systematics of the fits see [l3[ and 
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Lat 



ro 



a[fm] 



7.5 6.2875(10) 0.07952(1) 



10.0 8.6602(8) 0.05812(1) 



10.0 8.6602(8) 0.05812(1) 



12.5 10.916(3) 0.04580(1) 



R 



{Si} T T[fm] t s size N s N t #meas 



7-20 A 



9-27 A 



9-27 B 



11 - 29 A&B 



6 0.477 

10 0.795 

14 1.113 

18 1.431 



4 38 3 36000 



40 
42 3 
54 3 



3 



1500 4400 
3000 6468 
9000 11176 
18000 6512 



8 0.465 6 40 3 48000 3000 



10 0.581 4 50 3 



14 0.814 
18 1.046 



6 56 
4 54 



3000 
6000 
12000 



1272, 1272, 1296* 
2352, 2544, 2568 
6384, 6216, 6480 
8664, 8304, 8472 



6 0.349 4 48 3 48000 1500 

8 0.465 6 48 3 3000 

10 0.581 4 50 3 6000 

14 0.814 6 56 3 12000 



2000, 2000* 
2000, 6000 
2000, 7960 
2000, 2000 



8 0.366 6 48 3 36000 2000 
10 0.458 4 50 3 3000 



14 0.641 6 56 3 
18 0.824 4 72 



3 



6000 
12000 



1000 
4000 
7080 
2080 



TABLE II: Run parameters of the simulations. The number of measurements marked with * corresponds to the R values 
9 — 15, 17 — 21, 23 — 27 and the ones marked with ** to 9 — 15, 17 — 27 respectively. 



E. Finite volume effects 



below our statistical errors. 



Finite volume effects are possible due to finite spatio- 
temporal extents of the lattices. To control the effects 
due to finite temporal extents, we use four different ex- 
tents to make sure that any systematic error in our energy 
values is much lower than our statistical errors. 

The dominant corrections due to finite spatial ex- 
tents are around-the- world glueball exchanges of the form 
a(R) exp(— niQ{L — R)) , where ma is the glueball mass 
and L is the spatial extent of the lattice. These are rel- 
evant for large values of R. While glueball masses for 
3-dimensional SU(2) lattice gauge theory has been stud- 
ied in [23l l24j. the mixing coefficient a(R) is not known. 
To get an estimate for this coefficient, we did runs at dif- 
ferent lattice volumes at (3 — 5 (for which ma has been 
measured in [HI), varying L by a factor two. Since we 
did not see any finite size effects in the energies, we con- 
cluded, that in the range considered, a(R) is a number 
of order one. 

Using interpolation, we obtain from the data in 24], 
the glueball masses at (3 = 7.5, 10.0 and 12.5 to be 0.856, 
0.671 and 0.536 respectively in lattice units. Since our 
data at large R, with systematic errors under control is 
only at (3 = 7.5, we concentrate on that (3 value. The only 
lattice where this correction matters is the L = 42 lattice 
on which we measure Wilson loops with T = 14 and R 
between 7 — 20. On this lattice, the correction ranges 
between 10~ 9 and 6.7 x 10~ 9 for R between 18 - 20. 
The error bars for the wilson loop corresponding to the 
second excited state are of similar magnitude for R = 19 
and 20 while for the first excited state, R = 20 has a 
similar error bar. We therefore expect that there can 
be some correction due to the glueball for these energy 
values while for all other values, the corrections are well 



III. SIMULATION RESULTS 

In this section we discuss the results of our simulations, 
that are tabulated in table Mil We compare to the full 
NG spectrum, eq. (JTJ) and to leading order (LO) and next 
to leading order (NLO) models which are obtained by 
truncating the expansion of the square root of |T]) in 1/R 2 
at leading order and next to leading order respectively. 
The LW and PS type string theories give identical results 
to NLO order. The curves in the plots were drawn using 
the string tension at (3 = 12.5. All the results have been 
rescaled so that they are visible on a single plot. 

A. Energy states 

We use the groundstate to determine the string ten- 
sion and fix an additive constant Vq, appearing in the 
potential, by fitting to the form 



V(R) = aRJl 



12 a R 2 



Vn 



(7) 



The results of the fits are shown in table IIVI where also 
the values from [18| are shown for comparison. We see 
that all results agree well within error bars. Also a and 
Vq from both operator sets agree with each other. 

In figure [5] the results for the total energies are shown 
against R/tq, together with the predictions of the NG 
spectrum. The energies are rescaled such that = n. 
The results for the groundstate are completely in agree- 
ment with the NG predictions. In this case only the curve 
of the full spectrum is shown since deviations at different 
orders of the expansion are not visible on this scale. 
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T(18, 14) 


4T * 


1 (18, 14) 
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T(18, 14) 
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0.4251(1) 


0.42992(8) 


0.753(1) 


0.773(2) 


0.91(7) 




0.328(1) 


0.302(3) 


0.49(7) 
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0.4660(1) 


0.47226(9) 


0.768(1) 


0.784(2) 


0.96(6) 




0.302(1) 


0.284(2) 


0.50(6) 
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0.5064(2) 


0.5135(1) 


0.785(1) 


0.797(3) 


0.99(5) 




0.279(1) 


0.267(3) 


0.49(5) 
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0.5464(2) 


0.5464(2) 


0.806(1) 


0.804(2) 


1.01(4) 




0.259(1) 


0.257(3) 


0.46(4) 
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0.5863(2) 


0.5862(3) 


0.828(1) 


0.826(2) 


1.03(3) 




0.242(1) 


0.240(3) 


0.43(7) 
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0.6259(3) 


0.6259(3) 


0.853(1) 


0.850(2) 


1.04(2) 




0.227(1) 


0.224(3) 


0.41(2) 
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0.6655(3) 


0.6654(4) 


0.879(1) 


0.876(3) 


1.06(2) 




0.214(1) 


0.211(3) 


0.40(2) 
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0.7049(3) 


0.7047(5) 


0.907(1) 


0.904(3) 


1.08(2) 




0.202(1) 


0.199(3) 


0.38(2) 
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U. 1 442(4) 


U. (441(0) 


0.930(1) 


0.933(3) 


1.10(2) 




0.192(1) 


n i on/ o \ 
0.189(3) 


0.36(2) 
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0.7835(4) 


0.7833(6) 


0.966(2) 


0.963(3) 


1.13(1) 




0.183(2) 


0.180(4) 


0.34(1) 
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0.8228(5) 


0.8225(7) 


0.998(2) 


0.995(4) 


1.15(1) 




0.175(2) 


0.172(4) 


0.33(1) 




18 


0.8620(5) 


0.8617(7) 


1.030(2) 


1.027(4) 


1.18(1) 




0.168(2) 


0.165(4) 


0.32(1) 
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0.9012(5) 


0.9009(8) 


1.064(3) 


1.061(5) 


1.21(1) 




0.162(3) 


0.160(5) 
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0.9404(6) 


0.9401(8) 


1.097(4) 


1.095(6) 


1.24(2) 




0.157(4) 


0.155(6) 
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m /-to i a \ 

1 (18, 14) 


4T * 


T(18, 14) 


3T * 


T(10,8) 


3T* 


T(18, 14) 
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0.3160(1) 


0.3207(1) 


0.571(1) 


0.563(4) 


0.772(4) 


0.770(9) 


0.255(1) 


Best 


Best 
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0.3599(1) 


0.3596(5) 


0.587(1) 


0.567(6) 


0.770(3) 


0.757(7) 


0.227(2) 


estimate 


estimate 
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0.4032(2) 


0.4026(8) 


0.609(2) 


0.588(5) 


0.780(3) 


0.769(7) 


0.206(2) 


is 


is 
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0.4460(2) 


0.445(1) 


0.634(2) 


0.613(5) 


0.794(4) 


0.780(8) 


0.188(2) 


E\ — E» 


E2 — Eo 
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C\ A O O A f 0\ 

0.4884(3) 


0.487(1) 


O.oo3(3) 




0.809(2) 




0.175(3) 
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0.5309(3) 


0.529(2) 


0.694(3) 




0.834(3) 




0.163(3) 
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0.5733(4) 


0.571(2) 


0.726(4) 




0.860(3) 




0.152(4) 
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0.6153(5) 


0.613(4) 


0.757(7) 




0.889(4) 




0.142(7) 
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t~\ r 1 r T r / o\ 

0.6575(8) 


0.655(4) 


0.79(1) 




0.919(6) 




0.14(1) 
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0.699(2) 


0.697(4) 


0.84(3) 




0.954(10) 




0.14(3) 
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T(14, 10) 
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0.2540(2) 


0.2535(5) 


0.472(1) 


0.444(9) 


0.612(3) 




0.218(1) 


Best 


Best 
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0.2953(3) 


0.2945(9) 


0.488(2) 


0.459(9) 


0.613(3) 




0.192(2) 


estimate 


estimate 
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0.3359(4) 


0.335(1) 


0.510(2) 


0.481(8) 


0.625(3) 




0.173(2) 


is 


is 
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0.3761(6) 


0.374(2) 


0.536(3) 


0.512(12) 


0.644(4) 




0.158(3) 


E\ — Eo 


E2 — Eq 
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0.4162(9) 


0.414(4) 


0.566(5) 




0.663(5) 




0.147(5) 








26 


0.456(1) 


0.453(6) 


0.593(9) 




0.688(7) 




0.134(9) 








29 


0.496(2) 


0.492(10) 


0.622(16) 




0.707(11) 




0.122(16) 







TABLE III: Results of the simulations. The data columns used in the plots have been denoted by a * on their headers. The 
left column corresponds to the naive value of loops, obtained from Wilson loops with the temporal extents given in brackets 
in the header. The right column corresponds to the value from a fit to the form (JJ),©, wherever this was possible and the 
resulting values are expected to have negligible systematic errors. The header 4T means, that the values were obtained using 
fits to all possible combinations of the four different temporal extents, and the header 3T means that combinations from the 
three lowest temporal extents were used. Where these fits did not work and the data may contain significant systematic errors, 
we have only indicated what our best estimate is. The naive data serves both for future referencing as well as estimating the 
corrections due to finite temporal extents. They give useful upper bounds. 



For the excited states we were only able to obtain cor- 
rected results over a limited range of R, except for E\ at 
(3 = 7.5 where we have results over the full range. 



We see that the corrected values for /3 = 7.5 follow 
the Nambu curves quite well. In the region where we are 
able to obtain corrected energies for all values of (3, we 
see good agreement between the different j3 values. At 
smaller values of R where the difference between the NG, 
LO and NLO curves are visible, the data seems to favor 
the NG curve. 



B. Energy differences 

We now turn to the energy differences which are more 
sensitive to subleading properties of the flux tube. In ad- 
dition they have the advantage that the constant Vb does 
not contribute. The results for the energy differences 
AE10 and AE20 are shown in figure 

We were able to obtain corrected values using eq.© 
only for AE\q at (3 = 7.5. That data set seems to follow 
the NG curves quite nicely. For AE20 , our best estimates 
are the naive values from eq.© with T a — 14 and Tb = 
18. Nevertheless we expect very little systematic effects 
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FIG. 5: Results for the total energies E n . The 0' s are the 
values for f3 = 7.5, D's for (3 = 10.0 and O's for f3 = 12.5. V 
are corrected values for /3 = 10.0 and A for f3 = 12.5. The 
results have been rescaled, such that E n ° — n. The labeling 
of the lines are as defined in the text. 



in these results as the physical temporal extents for these 
loops are > lfm. Again we see that the NG curve is 
favored by the data. 

For AE W at = 10.0 and 12.5, our best estimates 
come from the difference E\ — Eq with the errors being 
calculated in quadrature. Even though this gives larger 
error bars, especially at (3 = 12.5, we see that both data 
sets are consistent with (3 — 7.5 values. 

For AE20 at (3 = 10.0 and 12.5, the physical temporal 
extents of the loops are not large enough that the higher 
order corrections, unaccounted for in equations (|3]) and 
(O, are negligible. Nevertheless, at (3 — 10.0, where it 
was possible at least partly to take into account the cor- 
rections, we see the trend of the data is to follow the 
NG curve. For j3 = 12.5, where even that was not possi- 
ble, the data contains systematic errors and we have not 
plotted it. 



AE20 




R/r 



FIG. 6: Top: Results for the energy difference AEio- Bot- 
tom: Results for the difference A.E20. The mapping of the 
curves and the points is the same as in the plot for the total 
energies. 



IV. CONCLUSIONS 

In this Letter we have looked at a variant of the multi- 
level error reduction scheme suitable for studying the ex- 
cited states of the QCD flux tube. Using this scheme we 
have looked at the excited states of the flux tube at three 
different couplings for qq separations between 0.5 and 1.7 
fm. 

Compared to measurements with the older method 
[l4l |. we have seen very significant improvements due to 
the new measurement process. It has been possible to 
increase the range of R by more than 50%. Since the er- 
ror in AE grows roughly proportional to R 2 , this would 
be very hard to achieve by increasing computing power 
alone. Compared to [l4j], the errors on E\ and AE^o 
have been reduced by a factor between 2 and 3 in the 
overlapping range of R. This would also be very hard 
to achieve by increasing statistics alone. Both of these 
observations clearly demonstrate the superiority of the 
new measurement process over the older one. Moreover 
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using the older method we did not have a signal for the 
second excited state whereas we do have a reasonable 
measurement now at least for one (3 value. 

Our results seem to indicate that in this scheme, with 
our basis, upto the values of R we have considered, one 
needs temporal extents of about 1 fm to make sure that 
the corrections due to the higher states are well under 
control. Wherever this criterion has been met, we have 
seen that the data follows the NG curve. 

The corrected data sets for different (3 values seem to 
fall on top of each other indicating that there is very little 
effect due to finite lattice spacing. As such it seems to be 
much more important to go to larger temporal extents 
than finer lattices. 

Unfortunately our data is still not good enough to dis- 
tinguish between NG, LW and PS type string theories as 
that would require a sensitivity at the level of R~ 5 . De- 
viations from the NG predictions at higher orders have 
recently been reported in [25[ for gauge duals of random 
percolation problems. 

This was our first attempt to show how to combine 
the two powerful techniques of multi-level error reduction 



and sophisticated sources for Wilson loops. Using this we 
have been able to go to much larger values of R and T 
than was possible before. For ruling out different string 
models, the error bars will have to be reduced by at least 
an order of magnitude. That at the moment looks like a 
project for the future. 
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p 


7.5 set A [18] 


10.0 set A 10.0 set B [18] 


12.5 set A 12.5 set B [181 


a 
V 


0.03867(7) 0.038566(6) 
0.1730(5) 


0.0206(2) 0.0210(4) 0.020606(4) 
0.145(3) 0.139(6) 


0.0128(4) 0.0129(4) 0.012742(17) 
0.124(4) 0.124(4) 



TABLE IV: Results for the string tension a and the regularization constant Vb- The values [l8| ] are the reference results for 
the string tension. 



